VLA Detection of the Ionized Stellar Winds Arising from Massive 
Stars in the Galactic Center Arches Cluster 
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ABSTRACT 

The Galactic center Arches stellar cluster, detected and studied until now only in the near- 
infrared, is comprised of at least one hundred massive (M* > 20 M Q ) stars. Here we report 
the detection at centimeter wavelengths of radio continuum emission from eight radio sources 
associated with the cluster. Seven of these radio sources have rising spectral indices between 4.9 
and 8.5 GHz and coincide spatially with the brightest stars in the cluster, as determined from 
JHK photometry and Bra and Br7 spectroscopy. Our results confirm the presence of powerful 
ionized winds in these stars. The eighth radio source has a nonthermal spectrum and its nature 
is yet unclear, but it could be associated with a lower mass young star in the cluster. 

Subject headings: Galaxy: center — stars: winds — stars: mass loss 
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1. Introduction 

High-resolution observations of the Galactic 
center at near-infrared wavelengths first revealed 
the remarkable Arches stellar cluster located at 
i=(F.l2, b=0?02, -11' (or 27.5 pc) in projection 
from Sgr A* (Nagata et al. 1995; Cotera et al. 
1996). Here, we assume a distance to the Galac- 
tic center of 8.0 kpc (Reid 1993). Nagata et al. 
(1995) identified the 14 brightest stars in this 
cluster with JHK photometry and Br7 and Bra 
hydrogen recombination lines, showing that these 
stars had the characteristic colors and emission 
lines of Wolf-Rayet (WR) and He I emission-line 
stars. The strength of the Bra and Br7 emis- 
sion lines suggests that these stars are losing mass 
at rates of M~2xl0~ 5 M yr" 1 . Using K-band 
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spectroscopy, Cotera et al. (1996) also found 12 
stars in the cluster with spectra consistent with 
late-type WN/Of types, and which have M~l-20 
x 10~ 5 M Q yr" 1 , and Voo -800-1200 km s _1 . 

Recent observations of the Arches cluster in- 
dicate that the number of massive stars is much 
greater than suggested by the initial studies. Ser- 
abyn et al. (1998) estimate that at least 100 clus- 
ter members have masses > 20 M , and Figer et 
al. (1999) calculate that the total cluster mass ex- 
ceeds 10 4 M Q with as many as 160 O-stars and a 
stellar density of p ^3xl0 5 M Q pc~ 3 . These pa- 
rameters make the Arches cluster one of the dens- 
est and most massive young stellar clusters in the 
Galaxy. In addition, the ionizing flux (estimated 
to be a few x 10 51 s _1 ) is sufficient to ionize the 
surrounding Arched Filaments H II complex (Ser- 
abyn et al. 1998; Lang et al. 2001). The relative 
placement of the Arches cluster and the filamen- 
tary nebular features is depicted in Figure 1. 

Thermal radio emission from ionized stellar 
winds is detectable at radio wavelengths and arises 
from the outer parts of the wind envelope. The 
classic theory of Panagia & Felli (1975) and Wright 
& Barlow (1975) predicts that the spectrum of ra- 
dio emission from a stellar wind is oc jy+°- e for 
a spherically symmetric, isothermal, stationary 
wind expanding at a constant velocity. Previous 
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surveys made with the Very Large Array (VLA) 
of the National Radio Astronomy Observatory 4 
have detected radio emission from OB supergiant 
and Wolf-Rayet stellar winds (Abbott et al. 1986; 
Bieging et al. 1989). Recently, Lang et al. (1999) 
detected six radio point sources in the vicinity of 
a similar Galactic center cluster, the Quintuplet, 
which have spectra consistent with the theoretical 
predictions and close positional correspondances 
with known near-infrared stellar sources. Here, 
we present multi-frequency VLA observations of 
the Arches cluster in an attempt to detect stellar 
wind emission from the individual stars known to 
be members of this dense cluster. 

2. VLA Continuum Observations 

VLA multi-frequency continuum observations 
were made of the Arches cluster in July 1999 (A- 
array) at 4.9 GHz and 8.5 GHz (4 hours each), 
and in October 1999 (BnA array) at 43.3 GHz (3 
hours). The phase center of all observations was 
the position given for the bright, central source in 
Nagataetal. (1995): a,<5 (J2 ooo) = 17 45 50.41, -28 
49 21.8. Standard procedures for calibration, edit- 
ing, and imaging were carried out using the AIPS 
software of NRAO. At all wavelengths, 3C286 was 
used for flux density calibration, and at 8.5 and 
4.8 GHz, 1751-253 was used as a phase calibra- 
tor. At 43 GHz, fast switching between the bright 
source Sgr A* and the phase center of the obser- 
vations (an angular distance of 12'5) was carried 
out on a 90 second cycle in order to closely track 
the phase variations. The resulting spatial resolu- 
tions and rms noise levels in the images are 0'.'79 
x 0'.'33, PA=12?8 and 45 ^Jy beam" 1 at 4.9 GHz; 
0741 x 0716, PA=2?8 and 25 ^Jy beam" 1 at 8.5 
GHz; and 0729 x 0707, PA=49?2 and 300 
beam" 1 at 43.3 GHz. 

3. Results 

Figure 2 shows the 8.5 GHz image of the ra- 
dio sources within a 16" (0.7 pc) region centered 
on the bright source at the cluster center. The 
crosses represent the positions of 9 of the 14 stars 
in Nagata et al. (1995) which showed Br7 and Bra 

4 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. 




Fig. 1. — Location of Arches stellar cluster and 
the surrounding Arched Filament H II nebula. 
The contours represent VLA 8.5 GHz continuum 
emission at levels of 7.5, 15, 20, 25, 37.5, 50, 75, 
87.5 mJy beam -1 , with a resolution of 778 x 676, 
PA=-0?7, from Lang et al. (2001). 

recombination lines. These positions have been 
shifted by 0^05 (0766) to correct for what appears 
to be a systematic offset to the West. At 8.5 GHz, 
eight sources are detected above a 6a level and are 
referred to as AR1— 8. Seven of the sources (with 
the exception of AR6) appear to be associated 
with the stellar positions plotted in Figure 2. AR7 
also shows a larger displacement (Aa=075) from 
the corresponding stellar position than the rest of 
the radio sources. However, the absolute uncer- 
tainty in the near-infrared positions is ~2" or less, 
so an association is quite possible within these lim- 
its. Since the a priori probability of finding a 8.5 
GHz background source with flux density greater 
than 0.1 mJy in the region considered is ~0.005, 
we conclude that all the radio sources detected are 
indeed associated with the Galactic center Arches 
cluster. 

The positions and properties of AR1— 8 are 
listed in Table 1. Four of these sources are de- 
tected in the 4.9 GHz observations (AR1, AR4, 
AR6, AR7), while the central source AR1 is ad- 
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Fig. 2. — VLA 8.5 GHz continuum image of the 
inner 16" (0.7 pc) of the Arches cluster, with a res- 
olution of Of! 41 x 0'/16, PA=2?8. This image has 
been primary beam corrected and has an rms of 
25 /xJy beam -1 . The contour levels of 0.075, 0.15, 
0.25, 0.5, 1.0, 1.5, 1.75 mJy beam -1 correspond to 
3, 6, 10, 20, 40, 60, 70 a and the crosses corre- 
spond to the positions of the nine near-infrared 
sources of Nagata et al. (1995) detected in Br7 
and Bra lines. These positions are only known to 
an absolute accuracy of ^2" (or 0?15), and here 
have been shifted by 0"66 (or 0?05) to the West to 
align the sources more closely. 

ditionally detected in the 43.3 GHz observations. 
At 43.3 GHz, AR1 has a flux density of 3.1±0.1 
mJy, consistent with the 4.9/8.5 GHz spectral in- 
dex (a^+0.35). Spectral indices between 4.9/8.5 
GHz are calculated for the sources detected at 4.9 
GHz. Upper limits for the flux densities at 4.9 
GHz can be estimated for the remaining sources, 
with resulting lower limits for the spectral indices. 
Seven of the sources show rising spectra in the 
range of +0.3 to +0.9 with the exception of AR6, 
which has a nonthermal spectral index (a=— 0.7). 
The sources AR2— 8 are point sources with angu- 
lar sizes < 0'.'25. The deconvolved size of AR1 at 
8.5 GHz is O'.'ll (0.005 pc), implying a brightness 
temperature of ^3000 K, a reasonable value for an 
ionized wind. 



4. Discussion 

4.1. The Nature of the Radio Sources 

The spectral indices of the detected sources 
are consistent with values expected from stellar 
wind emission. Several of the sources have spec- 
tra which are flatter (a ^+0.3) than the predicted 
a=+0.6, and AR6 has a negative spectral index 
(a ~— 0.7). However, these spectra do not rule 
out stellar wind detections. About 25% of stel- 
lar winds arising from massive stars have been 
found to have flat or nonthermal spectra in the ra- 
dio regime (Bieging, Abbott & Churchwell 1989). 
The precise mechanism that produces these non- 
thermal spectra is not established, but internal 
shocks in the wind of a single star (White 1985) 
or the presence of interacting winds in close bi- 
nary systems have been proposed (Contreras & 
Rodriguez 1999). For a typical late-type WN star 
at the Galactic center, the expected flux density 
at 8.5 GHz arising from a stellar wind is ^250 /iJy 
based on the assumed values of Vqo^IOOO km s -1 , 
M~5xl0 -5 M Q yr -1 and T-10 4 K. This value is 
similar to the flux densities we measure at 8.5 GHz 
for AR2— 8. Figure 2 also shows an excellent po- 
sitional correspondence between the radio sources 
and the stellar sources detected by Nagata et al. 
(1995). Table 2 lists the spectral types of these 
stars as determined by Cotera et al. (1996). The 
stars have been classified as late-type WN and Of 
stars, all of which are expected to be losing mass 
at high rates from their surfaces. As noted above, 
such stars often exhibit flattened or nonthermal 
spectral indices due to the presence of a nonther- 
mal component. Therefore, the range of spectra 
we detect is consistent with predictions for late- 
type WN stars. 

The radio source AR6 is displaced by Aa=l" from 
the corresponding stellar position of Nagata et al. 
(1995). At 8.5 GHz, there is a 75 /xJy source 
(3cr) closer to the stellar position, which may be 
the correct identifcation. The stellar source at 
this position has been classified as a WN8 star 
by Cotera et al. (1996). If this displacement is 
real, we have then to consider the possibility that 
AR6 is associated with another star in the cluster. 
T Tauri stars are sometimes found to have asso- 
ciated nonthermal emission, usually interpreted 
as gyrosynchrotron produced by electrons acceler- 
ated in situ by magnetic reconnection flares near 
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the stellar surface (Andre et al. 1988; Rodriguez, 
Anglada, and Curiel 1999). If this is the case, we 
may be detecting for the first time a young, low 
mass member of the cluster. 

4.2. Mass Loss Rates 

The flux densities at 8.5 GHz were used to cal- 
culate a mass-loss rate based on Panagia & Felli 
(1975): 



1O- 5 M yr~ 



A - o.52 (^l) 3/4 ( HSS ) (J- 

\mjy J V 10 3 km s _1 / \kpc 

(1) 

where Sg.s is the flux density of the source at 8.5 
GHz, Vqo is the terminal velocity of the wind, and 
d is the distance to the source (8 kpc). We have 
assumed an electron temperature of T*=10 4 K, 
Z=l, and a mean molecular weight, /i=2, due to 
the enrichment in heavy elements of the late-type 
WN stars (Lcitherer et al. 1997). The mass loss 
rate of AR1 has the largest value, M=1.7 x 1CT 4 
M Q yr _1 , comparable only to stars with extreme 
values for M, such as in the R136 stellar cluster at 
the center of 30 Doradus (de Koter et al. 1997). 
The values of M for AR2— 8 are in the range of 
3—4.5 x 10~ 5 M Q yr _1 , consistent with values 
found for for late type WN stars by Leitherer et 
al. (1997): M~2. 5-3.9 xlCT 5 M yr" 1 . 

4.3. The Arches Cluster and the Interstel- 
lar Medium 

The rich population of massive stars in the Ra- 
dio Arc region, which includes both the Arches 
and Quintuplet clusters (Figer et al. 1999), sug- 
gests that the interstellar environment should be 
strongly influenced by the energetic processes re- 
lated to such massive star formation and evo- 
lution. Ionization of the molecular material in 
this region and formation of the unusually shaped 
Arched Filaments and Sickle H II regions can be 
accounted for by the ionizing fluxes of each clus- 
ter (Lang et al. 1997, 2001). In addition, inter- 
actions between the winds in the densely-packed 
Arches cluster may give rise to a shock-heated dif- 
fuse "cluster wind" , which could have a tempera- 
ture as high as 10 7 K (Canto, Raga & Rodriguez 
2000). This emission is predicted to be detectable 
in the X-ray regime with the Chandra X-ray Ob- 
servatory. Studies of the interaction between the 



massive stars in the Arches and the surrounding 
interstellar environoment are currently underway 
at a number of complementary wavelengths. 

5. Conclusions 

We present the following results from our 4.9, 
8.5, and 43.3 GHz observations of the Arches Clus- 
ter: 

(1) Eight radio sources are detected at 8.5 GHz, 
four sources at 4.9 GHz, and one of them is de- 

3 / 2 tected also at 43.3 GHz. Seven of the eight de- 
tections show rising spectral indices between 4.9 
and 8.5 GHz in the range +0.3 to +0.6, while one 
source shows a nonthermal spectrum (—0.7). All 
of the detections are consistent with stellar wind 
emission from mass-losing stars. However, the 
nonthermal source may be tracing a young, low 
mass member of the cluster. This is the first de- 
tection of emission associated with the members 
of the Arches cluster at a band other than the 
near-infrared. 

(2) Seven of the eight radio sources are coinci- 
dent in position with the Br7 and Bra emission- 
line stellar sources of Nagata et al. (1995), which 
correspond to the late WN and Of stellar types 
as classified by Cotera et al. (1996). These po- 
sitional correspondences further confirm that the 
ionized winds from stars in the Arches cluster arc 
detected. 

(3) Based on the radio flux density, AR1 has a 
mass-loss rate of ^2xl0~ 4 M Q yr _1 , comparable 
only to that of the stars in 30 Doradus (de Koter 
et al. 1997). Mass-loss rates for AR2— 8 are found 
to be in the range 3— 4.5xl0 -5 M yr _1 , similar 
to the large mass loss rates for other late-type WN 
systems (Leitherer et al. 1997). 
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Table 1 

Parameters of the Radio Sources (AR1-8) in the Arches Cluster 



Source RA (J2000) t DEC (J2000)* Flux Density (mJy) Spectral M 



Name 


(h m s) 


(° ' ") 


^=8.5 GHz 


^=4.9 GHz 


Index 


(M Q yr- 


- 1 ) 


AR1 


17 45 50.42 


-28 49 22.0 


1.70±0.05 


1.40±0.03 


+0.35±0.04 


1.7x10" 


-4 


AR2 


17 45 50.39 


-28 49 21.0 


0.23±0.02 


<0.13 


>+0.9 


3.9x10- 


-5 


AR3 


17 45 50.20 


-28 49 22.0 


0.17±0.02 


<0.13 


>+0.4 


3.2x10" 


-5 


AR4 


17 45 50.47 


-28 49 19.2 


0.23±0.02 


0.16±0.03 


+0.65±0.13 


3.9x10- 


-5 


AR5 


17 45 50.57 


-28 49 17.2 


0.16±0.02 


<0.13 


>0.3 


3.0x10- 


-5 


AR6 


17 45 49.76 


-28 49 25.7 


0.27±0.02 


0.40±0.03 


-0.71±0.08 


4.5x10- 


-5 


AR7 


17 45 50.83 


-28 49 26.1 


0.25±0.02 


0.21±0.03 


+0.31±0.05 


4.2x10- 


-5 


AR8 


17 45 50.45 


-28 49 31.7 


0.20±0.02 


<0.13 


>+0.7 


3.6x10" 


-5 



^The positional accuracy in RA is ±0^01 
^The positional accuracy in DEC is ±0'.'l 



Table 2 

Near-Infrared Counterparts of Arches Radio Sources 



Radio 


Nagata et al. 


Cotera et al. 


Spectral 


Source 


Source # 


Source # 


Type* 


AR1 


8 


8 


Of/WN9 


AR2 


7 


6 


WN8/9 


AR3 


4 


9 


Of/WN9 


AR4 


10 


5 


Of/WN9 


AR5 


11 


2 


WN8 


AR6 


It 


13 


WN8 


AR7 


14 


11 


Of/WN9 


AR8 


9 







*From Cotera et al. (1996). 

' Identification with near-infrared source not firm. 
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